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Abstract 


Within the MSC/NASTRAN DMAP module TRD1, solving physical (coupled) or modal (uncoupled) transient 
equations of motion is performed using the Newmark-Beta or mode superposition algorithms, respectively. For 
equations of motion with initial conditions, only the Newmaric-Beta integration routine has been available in 
MSC/NASTRAN solution sequences for solving physical systems and in custom DMAP sequences or Alters for 
solving modal systems. In some cases, one difficulty with using the Newmark-Beta method is that the process of 
selecting suitable integration time steps for obtaining acceptable results is lengthy. In addition, when very small step 
sizes are required, a large amount of time can be spent integrating the equations of motion. For certain aerospace 
applications, a significant time savings can be realized when the equations of motion are solved using an exact 
integration routine instead of the Newmaric-Beta numerical algorithm. In order to solve modal equations of motion 
with initial conditions and take advantage of efficiencies gained when using uncoupled solution algorithms (like that 
within TRD1), an exact mode superposition method using MSC/NASTRAN DMAP has been developed and 
successfully implemented as an enhancement to an existing coupled loads methodology at the NASA Lewis Research 
Center. 
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Introduction 


Linear transient response analyses can be performed with MSC/NASTRAN in either physical or modal space. In 
many aerospace applications, the size of the physical system, in terms of the number of DOF, is large. Hence, the 
cost associated with directly solving for transient responses is prohibitive. Using the mode-superposition method, the 
physical system is reduced to a modal system. The number of modal DOF required for accurate transient analysis 
is usually much less than the number of physical DOF, and the resulting equations of motion are uncoupled. Hence, 
solving for the system response in modal space is usually much more efficient. 

For most aerospace applications, the system is free-free and is assumed to reach steady-state equilibrium before the 
transient analysis begins. Hence, initial conditions exist. Such is the case for STS liftoff and landing transient 
analyses. These analyses can be performed directly in physical space, but doing so is usually very costly due to the 
number of physical DOF. It is more efficient to solve for the system response in modal space. When using 
MSC/NASTRAN solution sequences to solve modal equations of motion without initial conditions, the uncoupled 
solution algorithm [1] within DMAP module TRD1 is used. Unfortunately, solving modal equations of motion with 
initial conditions using MSC/NASTRAN solution sequences is not possible. Solving for modal transient responses 
when initial conditions are present can be done using custom DMAP or altered MSC/NASTRAN solution sequences. 
The NASA Lewis Research Center has had the capability for solving such systems for many years via custom DMAP 
sequences [2], The only limitation is that because initial conditions are present, the logic within module TRD1 routes 
execution to the coupled solution algorithm [1]. This algorithm is a variation of the Newmark-Beta method [3], The 
Newmark-Beta method is an unconditionally stable algorithm [4]. However, the algorithm has been observed in some 
cases, as with STS liftoff analyses, to be computationally inefficient. While the routine will converge to an answer, 
very fine integration intervals are usually required to obtain accurate results. This translates into lengthy analysis 
times and high costs. 


In order to take advantage of the uncoupled solver within module TRD1, a method was developed whereby the total 
system modal response is solved for as the sum of two modal responses. Due to the linearity of the system, the input 
loads can be represented as the sum of two applied loads. One load is the initial load applied to the system at time 
t=0.0. The second load is the difference between the total and initial applied loads. This second (shifted) load can 
be regarded as the input load time history with the horizontal axis shifted by the value of the initial load. Thus, the 
shifted transient load at time t=0.0 is equal to zero. The modal responses due to the initial load are the steady-state 
elastic and rigid-body modal initial conditions of the system. The responses are solved for using custom DMAP. The 
modal responses due to the shifted transient loads are solved for using die uncoupled algorithm of module TRD1. 
This is possible since the zero load at time t=0.0 translates into zero initial conditions. Hence, the logic within TRD1 
allows for the use of the uncoupled solver. Since the responses are assumed linear, the total modal response is the 
sum of the shifted transient and steady-state responses. 

The theory detailing the new exact mode superposition transient response methodology is presented in the next 
section. Following that, the implementation of the solution algorithm within the existing NASA LeRC coupled loads 
methodology is explained. Lastly, the solution process is illustrated with a numerical example for which an analytical 
solution is obtained. 


Theory 

For a given aerospace system, the coupled loads analysis process begins by assembling component models to form 
a system model. Generally, damping is ignored at the component level, so die assembled physical equations of 
motion for the system over all time instants are 


[MJ[u a ] + [K aa ][u a ] - [P a ] (1) 

where the number of columns of the physical response and load matrices is equal to the number of integration time 
instants. The ith column of each response matrix corresponds to a solution vector at the ith instant of time. For many 
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applications such as STS analyses, the analysis is performed on a free-free system assumed to reach steady-state 
equilibrium. By performing a system level eigenvalue analysis, a set of h mass normalized system mode shapes, 
[<J> ah ], m obtained. These consist of rigid-body and elastic mode shapes. The number of associated system modal 
DOF are usually much less than the original number of physical DOF. Using the mode-superposition method, the 
physical DOF are expressed in terms of the modal DOF, or 

[u a ] - M> ah ]E h ] (2) 

and 

[u fl ] - M>ah]tf h ] (3) 


where 5 are the system modal DOF comprised of rigid-body and elastic partitions. 

Substituting the expressions of Eq. (2) and Eq. (3) into Eq. (1) and premultiplying the resulting equation by the 
transpose of the system mode shapes, the system equations of motion are transformed and reduced from physical 
to modal space. The resulting system modal equations of motion are 


[MhhlKhl + [KhJRhl - [P h ] 

where 


IMhJ - 

(5) 

[Khh] ' NahftK JWahl 

(6) 

PU " Wal/tP.] 

(7) 


Note that [Mj^] [K hh ] are diagonal matrices. Since the system mode shapes are mass normalized, [M*^] is an 
identity matrix, and the terms of [K hh ] are the system eigenvalues. In most aerospace applications, the effects of 
damping are represented using proportional damping at the system modal level. Hence, the system damping matrix, 
[B^], is also diagonal. Adding the effects of system damping to Eq. (4) results in the final system modal equations 
of motion 


IMmPRJ + [B hh ]K h ] + [K hh ][? h ] * [P h ] (8 > 

Equation (8) is a set of h uncoupled single DOF system equations of motion. 

The modal input loads over all integration times are [P h ]. For the ith integration time instant, the modal load is {P^}. 
This load is equal to the steady-state equilibrium load at time t=0.0, {Pi}, plus the complementary (shifted) load, 
K%OT 
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{P h ‘} - {P h '} - {Ph ,sh } 


( 9 ) 


Hence, the modal loads over all time instants are 

[p h i - [pj] ♦ [p h sh i ( 10 > 

where each column of [P^] is {P^} , and the ith column of [P^* 1 ] is {Ph ,5h } • 11 is important to note that the first column 
of [PjJ h ] is {0}. 

Since the applied loads can be expressed as the sum of two loads as shown in Eq. (10), the system modal equations 
of motion are written as 


[MwJKhl ♦ \BM + [KhhlKhl - [P*] + [P h ‘] 


( 11 ) 


Linearity allows the total modal response to be equal to the sum of the responses due to each loading condition. The 
modal response due to [P£ h ] is a shifted transient response, and the modal response due to [P£| is a steady-state 
response. The two system modal equations of motion from which the two responses are obtained are 


IMhhlKf] + [B hh ][?’ h ] + tK hh ][lf] - IP*] 


( 12 ) 


and 


[M hh ][^‘] + [B hh ]R h “] + [K^lie**] - [Rjl 


( 13 ) 


Hence, the total modal displacements, velocities, and accelerations at each time instant tj are 


(4> - rif) + O 

(14) 


(15) 

- (tT) + > 

(16) 


Due to the fact that the initial loads in Eq. (12) are equal to zero, initial conditions are equal to zero for the shifted 
transient analysis. Hence, the shifted modal responses can be obtained using an exact uncoupled integration routine 
like that found within DMAP module TRD1. 

The initial conditions for the system are considered with the steady-state response of Eq. (13). The modal 
displacements (and velocities and accelerations) can be separated into r rigid-body and e elastic responses, or 
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( 17 ) 



Accordingly, Eq. (13) can be partitioned into rigid-body and elastic modal equations of motion, or 


[M „][«“] + [B„]Kr S ] + [KJlO = t P r‘] 


(18) 


and 


[M ee ][5 S e S ] + [B ee ][5 S e S ] - (K ee ][Cl = [ p e‘] 


(19) 


The elastic portion of the steady-state response is given by Eq. (19). Before the analysis begins, the system is 
assumed to reach steady-state equilibrium due to the initial loads. Hence, the elastic modal accelerations and 
velocities for all times are equal to zero in the steady-state condition: 


He] = Kf] = [0] 


( 20 ) 


Considering this, Eq. (19) is rewritten as 


MO - [ p e 'i 

The steady-state elastic modal displacements are solved for as 

iCl = l K ce]"'l P e l ] 


(21) 


( 22 ) 


Each column of steady-state elastic modal displacements is defined as |q 0 |. The values within |q,,| are constant over 
all time instants. 


The steady-state rigid-body modal response is found via Eq. (18). Since these DOF are associated with the rigid-body 
frequencies of the system. 


[B it ] = [KJ = [0] 


(23) 


and Eq. (18) is rewritten as 


[Mjir;*] = [p r ‘] 


(24) 
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Due to the normalization of the system mode shapes, [M rr ] is an identity matrix, and Eq. (24) reduces to 


£?] = [P r '] 


(25) 


The steady-state rigid-body modal accelerations are simply the accelerations due to the initial loads on the system, 
or 


Kr*] ■ [a] 


(26) 


where each column of [a] corresponds to the rigid-body modal accelerations, (a), due to the initial loads at time 
t-O.O. The values within |a) are constant over all times. The steady-state rigid-body modal velocities and 
displacements are found using the kinematic equations. The steady-state rigid-body modal velocities at each time 
instant tj are 



tila) 


(27) 


and the corresponding displacements are 

l£“l = ytflal (28) 

Substituting the steady-state responses into Eqs. (14), (15), and (16), the total modal responses at each time instant 
tj are: 


16 - 


O 

i£*i 


(q«l 


l«l» 


i*i sh i 


tjla) 


k 

itfi 

+ 4 

101 


\ll 


itfi 

► 

(a| 



► 

l°l 


(29) 


(30) 


(31) 


The exact mode superposition transient solution methodology presented above has been implemented within the 
NASA LeRC coupled loads methodology. 
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Implementation 


In order to solve for the total modal response of a system with initial conditions using MSC/NASTRAN uncoupled 
integration, a DMAP Alter can be written for a MSC/NASTRAN solution sequence, or a custom DMAP sequence 
can be developed. To implement the solution procedure within the NASA LeRC coupled loads methodology, the only 
modifications required were made to the custom DMAP sequence QTRAN [2]. The sequence is used to solve for 
system modal transient responses. As mentioned previously, the limitation of the DMAP had been that solving for 
modal responses of systems with initial conditions was only possible via the coupled integration routine within 
module TRD1. After implementing the theory presented in the previous section, solving for the responses of such 
systems is possible via either the coupled or uncoupled integrators available within TRD1. 

The exact mode superposition algorithm within the DMAP sequence QTRAN is as follows. To form the system 
modal mass, damping, and stiffness matrices as shown in Eq. (8), module GRAM is called. The physical loads 
applied to the system are assembled in a separate run and stored as matrix [PDT]. To form the shifted loads matrix, 
the first col umn of [PDT] is subtracted from all columns of [PDT]. Hence, all physical loads are shifted by the loads 
at time t=0.0, and the shifted transient load at time t=0.0 is equal to zero. Using the system mode shapes, the shifted 
and initial load matrices are transformed from physical to modal space. With the system modal mass, damping, 
stiffness, and shifted loads defined, module TRD1 is called to solve for the shifted modal response shown by 
Eq. (12). It is important to note that since the system modal equations of motion are uncoupled with zero valued 
initial conditions, the uncoupled integration routine is used. 

To solve for the steady-state system modal response defined by Eq. (13), several matrices are first generated. As 
described in the previous section, the system initial conditions are computed such that steady-state equilibrium is 
reached due to the initial applied loads at time t=0.0. The steady-state modal response is thus the system response 
due to constant loads applied to the system which are equal to the initial loads. To generate the steady-state response, 
two unknowns must be determined: the elastic modal initial conditions, {qj, and the steady-state rigid-body modal 
accelerations, {a}. The elastic modal initial conditions are found via Eq. (22). The steady-state rigid-body modal 
accelerations are found via Eq. (25). 

The steady-state modal elastic velocities and accelerations are zero because the system is assumed to reach 
equilibrium before the analysis begins. The first step in generating the steady-state rigid-body modal velocities and 
displacements is forming a vector of integration output times. The output times are stored in table FOL during 
original processing of the applied loads in a preceding run. Through a series of DMAP calls, the FOL table is read, 
and the time values are stored in a vector. This vector of times is used to generate the velocities of Eq. (27). By 
squaring each term of the time vector, the resulting vector of squared times is used to generate the displacements 
defined in Eq. (28). 

Given the steady-state elastic and steady-state rigid-body modal response matrices, a series of MERGE calls are made 
to form one matrix of steady-state modal responses. This matrix is then added to the matrix of shifted modal transient 
responses (see Eqs. (29), (30), and (31)) to form the final matrix of total modal responses of the system, [UHVF]. 
The modal responses are used throughout the remaining coupled loads methodology to perform physical data 
recovery operations. 


Numerical Example 


In order to exercise the new DMAP to the fullest extent and compare the results to an analytical (closed-form) 
solution, the following numerical example was developed. The example is that of a three DOF free-free system 
consisting of three masses and two springs. The system is shown in Fig. 1. Each mass is assigned a unit value, and 
the ratio of each stiffness to a mass is 100. Associated with the three DOF system is one rigid-body system mode 
and two elastic system modes. The mass and stiffness values for the system are such that the system circular natural 
frequencies are 0.0, 10.0, and 17.32 rad/sec. All modal DOF were retained, and the first DOF corresponds to the 
system rigid-body mode. 
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Applied to the system is the physical transient load shown in Fig. 2. It is a cosine function with a circular frequency 
of 15.0 rad/sec. Note that the load is nonzero at time t=0.0, and the load can be decomposed into two loads as shown 
in Fig. 2. Given the nonzero load at time t=0.0, the system is assumed to reach steady-state equilibrium before the 
analysis begins. Based on the theory presented previously, the system modal initial conditions are as follows: zero 
value rigid-body modal displacement and velocity and elastic modal velocities and accelerations, and nonzero value 
rigid-body modal acceleration and elastic modal displacements. 

The system modal responses for the numerical example were solved for using the new exact mode superposition 
method. The results were then compared to the closed-form analytical solution. Time histories of the modal 
accelerations, velocities, and displacements are shown in Fig. 3 through Fig. 5, respectively. Comparisons between 
maxima and minima data for the two solutions are shown in Table 1. From the figures and table, it is clear that the 
solutions obtained using the exact mode superposition algorithm are in exact agreement with the analytical solutions. 

Another numerical example using a real-world engineering problem was analyzed to compare the performances of 
the exact mode superposition methodology versus a solution methodology using a coupled integration routine. The 
results of this study are presented in [5]. 


Conclusion 


A solution algorithm has been implemented using MSC/NASTRAN DMAP whereby system modal equations of 
motion with initial conditions can be solved via the uncoupled integration routine wi thin DMAP module TRD1. The 
basis for the algorithm is that the total modal response due to applied loads on a system can be solved for as the 
superposition of a shifted transient response and a steady-state response due to initial loads. The exact mode 
superposition methodology has been implemented as an enhancement to the NASA LeRC coupled loads 
methodology. It has been shown via a numerical example that the exact mode superposition method is very accurate 
for solving system modal equations of motion with initial conditions. 


References 

[1] MSC/NASTRAN Programmer’s Manual, Vol. IV, The MacNeal-Schwendler Corp., Los Angeles, CA, 1985 
pp. 4.163-1 - 4.163-18. 

[2] Flanigan, C.C.; and Manella, R.T.: Advanced Coupled Loads Analysis Using MSC/NASTRAN. Proceedings 
of the 1991 MSC World Users’ Conference, The MacNeal-Schwendler Corp., Los Angeles, CA, 1991, Paper 
No. 14. 

[3] Newmark, N.M.: A Method of Computation for Structural Dynamics. J. Eng. Mech. Div., Am Soc Civ 
Eng., vol. 85, 1959, pp. 67-49. 

[4] Bathe, KJ.: Fmite Element Procedures in Engineering Analysis. Prentice-Hall, 1982. 

[5] Barnett, A.R., et al.: Solving Modal Equations of Motion with Initial Conditions using MSC/NASTRAN 
DMAP Part 2: Coupled vs. Uncoupled Integration. To be published in the Proceedings of the 1993 MSC 
World Users’ Conference, The MacNeal-Schwendler Corp., 1993. 


9 



Table 1. Numerical Example Results 




Ratio of Numerical Value to Anatytical Value 

Modal Response 

Extrema Response and Time of Occurrence* 



Modal DOF 1 

Modal DOF 2 

Modal DOF 3 


Minimum 

1.000 

1.000 

0.999 



(1.000) 

(1.000) 

(1.000) 

Acceleration 

Maximum 

1.000 

1.000 

0.999 



(1.000) 

(1.000) 

(1.000) 


Minimum 

0.999 

1.000 

0.999 



(1.000) 

(1.000) 

(1.000) 

Velocity 

Maximum 

0.999 

1.000 

0.999 



(1.000) 

(1.000) 

(1.000) 


Minimum 

1.000 

0.997 

1.000 



(1.000) 

(1.000) 

(1.000) 

Displacement 

Maximum 

1.000 

1.000 

1.000 



(1.000) 

(1.000) 

(1.000) 


* Ratio of time of occurrence values shown in parentheses. 




































► P(t) 



Figure 1. — System for numerical example. 
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Figure 2. — Applied load for numerical example. 



Figure 3. — Modal accelerations. 
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